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ABSTRACT 

Issues regarding the experimental implementation of PDE- based controllers are discussed in 
this work. While the motivating application involves the reduction of vibration levels for a cir- 
cular plate through excitation of surface-mounted piezoceramic patches, the general techniques 
described here will extend to a variety of applications. The initial step is the development of 
a PDE model which accurately captures the physics of the underlying process. This model is 
then discretized to yield a vector- valued initial value problem. Optimal control theory is used to 
determine continuous-time voltages to the patches, and the approximations needed to facilitate 
discrete time implementation are addressed. Finally, experimental results demonstrating the 
control of both transient and steady state vibrations through these techniques are presented. 
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1 Introduction 


An increasingly popular method for controlling structural vibrations is through the use of piezo- 
ceramic patches bonded to or imbedded in the structure. These pat ches exhibit the piezoelectric 
property that inplane strains are generated in response to an applied voltage. Depending upon 
the geometry of patch placement, location with respect to the structure’s neutral surface, and 
the method of excitation, this provides a mechanism for generating both inplane forces and/or 
bending moments in the underlying structure. 

The advantage of using such pat ches as actuators in many applications is due to the fact that 
they are lightweight, space efficient, relatively inexpensive, and provide a means of obtaining 
structural control without significantly changing the passive structural dynamics (they do not 
mass load the structure in the manner of a shaker or proof mass actuator). Due to their 
ceramic nature, they can also be molded in a variety of shapes so as to fit the structure under 
consideration. Moreover, rigid body torques and spillover effects are minimized due to the fact 
that the}’ are fully self contained and distributed in nature. Finally, they also exhibit the inverse 
piezoelectric property and hence generate a voltage in response to strains in the material. Hence 
a single patch or patch pair can be used for either sensing or actuation. This contributes to 
their efficiency and application in “smart material” structures. 

In addition to their utility in purely structural applications, they are also finding increasing 
use in structural acoustic and fluid/structure applications. Again, their advantage lies in the fact 
that they provide an efficient means of controlling the structure without significantly altering 
its passive dynamics. Through consideration of the coupling between the structure and the 
adjacent media, this provides a means ol controlling acoustic sound pressure levels or adjacent 
flow dynamics. 

A great deal of research in the last several years has been directed toward questions re- 
garding the modeling of piezoceramic patch interactions with underlying structures (see [11] 
and the references therein), strategies for determination of optimal patch location and number, 
and control techniques which utilize the patches as sensors and actuators. Due to the steady 
state periodic nature of the dynamics in many structural and structural acoustic systems, a 
large number of the current control methods are based upon frequency response input/output 
analysis. For example, Fuller et al. have employed a feedforward filtered X version of an 
adaptive LMS algorithm to control a simply supported plate using piezoceramic actuators [18]. 
While such techniques have proven quite successful for controlling steady state vibrations, they 
do not have the capability for direct control of transient responses. Other successfully im- 
plemented methods employing piezoelectric actuators to actively control structural vibrations 
include feedthrough techniques [17] and velocity feedback techniques [1, 24]. In general, these 
methods are based upon modal techniques and are designed to control purely steady state 
responses. An exception to this are the experimental results reported in [24] in which tran- 
sient plate vibrations, generated by an impact hammer, were reduced using an analog velocity 
feedback circuit. 

Similar studies have demonstrated the experimental success of using surface-mounted piezo- 
ceramic patches to reduce structure-borne noise in structural acoustic systems [19, 20]. The 
emphasis in these studies was again on using frequency input/output analysis to control steady 
state dynamics. 
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An alternative approach to controlling structural vibrations and sound pressure levels in 
structural acoustic systems is through the use of PDE-based feedback control methods. Analysis 
and numerical studies demonstrating these techniques for structural applications can be found 
in [6, 7, 8, 12] with corresponding results for structural acoustic systems given in [4, 5, 10]. These 
techniques start with an infinite dimensional PDE model for the system under consideration. 
When developing such models, care should be taken to incorporate not only the contributions 
due to the piezoceramic patches but also dynamics due to inexact boundary conditions [13], 
coupling with adjacent acoustic fields [5], as well as any other physical phenomena which affect 
the dynamics of the structure. In this setting, mathematical issues such as model welbposedness 
and approximation issues concerning simulations, parameter estimation and control can be 
addressed. 

By approaching the problems in this manner, one can avoid the difficulties caused in purely 
modal methods by patch contributions, coupling between components, and inexact boundary 
conditions. Moreover, by combining the PDE model with appropriate time-dependent feedback 
control theory, one obtains a method which is equally applicable for controlling transient or 
steady state vibrations. 

In this paper, the experimental implementation of such a PDE-based control method is con- 
sidered. While the motivating application involves the control of vibration levels for a circular 
plate through the excitation of surface-mounted piezoceramic patches, the general techniques 
described here will extend to a variety of applications. Following a brief discussion regarding 
the model and a Fourier/Galerkin scheme used to discretize it, relevant feedback control theory 
is discussed. In the discussion of the continuous and discrete time control results, two cases are 
considered; namely, the control of plate vibrations in the absence of a primary input force and 
the control of a plate driven by a periodic exogenous force. Implementation issues such as the 
effects of phase shifts and delays due to hardware are discussed and the experimental setup is 
briefly described. Finally, experimental results demonstrating the transient and steady state 
control I'esults are presented. These demonstrate the effectiveness of the PDE-based controller 
for this system and indicate the potential of these control techniques for reducing transient and 
steady state dynamics in other structural and structural acoustic systems. 

Finally, we offer comments on the nature of our contributions here to the literature. This 
paper does not contain any new theoretical results; it reports on our successful use of PDE- 
based methods in experiments at NASA Langley Research Center. The methods (which are 
theoretically sound) for parameter estimation and feedback control are based on approximation 
theory developed (by us and many others in the PDE control community) during the past 
several decades. All necessary theorems on convergence of finite dimensional parameter esti- 
mates, gains, filters, controls, observers, tracking variables, et cetera, needed for the examples 
treated here either have appeared or will soon appear in the research literature. Many of these 
theoretical results on PDE-based control have been largely viewed (especially in the engineer- 
ing community) as non-implementable and, hence, as somewhat irrelevant to applied scientists 
and engineers. The present manuscript refutes this notion and validates the practicality and 
importance of much of the theoretical efforts over the past years on PDE-based control. Our 
group is the first, to our knowledge, to provide irrefutable evidence that one can start from 
basic physical laws, derive careful infinite dimensional distributed parameter or PDE control 
models, and successfully implement PDE control methods (with the approximations necessary 
to obtain finite dimensional controls) which behave as theory and simulations predict. 
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2 Circular Plate Model 


The structure under consideration is a thin circular plate with s sets of piezoceramic patches 
bonded to the plate either singly or in pairs as depicted in Figure 2.1. Throughout this discus- 
sion, the radius and thickness of the plate are denoted by a and h , respect ively. Moreover, the 
density, Young’s modulus, Poisson ratio, Kelvin-Voigt and air damping parameters are given 
by Pi E , cd and 7 while w and g are used to denote the transverse displacement and external 
force on the plate. Finally, the region occupied by the unstrained neutral surface of the plate 
is indicated by To. 

Equations of motion for the plate can be determined from both Newtonian (force and 
moment balancing) and Hamiltonian (energy formulation) principles and we summarize both 
approaches here. The presentation will be for a general Kirchhoff plate with potentially non- 
axisymmetric responses. 



Figure 2.1. Thin circular plate with surface-mounted piezoceramic patches. 


2.1 Strong Form of Plate Model 

Considering first the model which derives from Newtonian principles, we let M r , Me , M T o denote 
internal moments and (M r ) p< ,,(Me) pe denote external moments generated by the piezoceramic 
patches. As detailed in [9, 29], for a structure with patch pairs activated out-of-phase, the 




internal moments are given by 
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where D = denotes the plate stiffness. The external moments generated by the s patch 

pairs have the form 

(M r ) pe = (M e ) pe = -i2>C?Ui(t)xi(r,0) , \ 1 (r.O) = 

2 = 1 

where Ui(t) is the voltage into the i th patch and Kf is a. parameter which depends on the geom- 
etry, piezoceramic material properties and piezoelectric strain constants (see [11] for details). 
The internal and external moments can then be combined to yield the general plate moments 

M r = Mr — (M r ) pe 
Me = Me ~ (M$) pt 
M r e M r e • 


1 , (r, 0) G i th patch region 

0 , elsewhere 


For a clamped plate with initial displacement iOo(r, 0) and velocity wi(r, 0), force and mo- 
ment balancing yield the equations 

d 2 w dw d 2 M r Id Mr 1 dMe 

ph ~dF ~ dr 2 ~r dr + r dr 0 < 0 < 2ir 

2d 2 M r e 2 dM r e 1 d 2 M e _ ' 0<r<« 

r drdO r 2 dO r 2 d9 2 g[ ' r ' > 

( du, (2-D 

< w(t, a , 0 ) = a i ^) = 0 



w(O,r,0) = tv o (r,0) , -^(0, r, 0 ) = t«i(r, 0 ) . 

as the strong form of the plate model (see [9, 27] for details). From an applications perspective, 
the following observations can be made regarding this model. 

1. It is first noted that in many applications, it is nearly impossible to maintain the truly fixed 
(zero displacement and slope) boundary conditions specified in the model (2.1). For the 

experimental plate which we used, parameter estimation results indicated minimal energy 
loss through the boundary conditions and an adequate fit of the model (2.1) was obtained 
(see the identification results in [2, 3]). For plates in which the boundary clamping is less 

secure, a model for imperfectly clamped boundaries such as that presented in [13] should 
be used in order to obtain a. model fit which is adequate for control applications. 
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2. As discussed in [11], the plate parameters p, D,// and cp are discontinuous due to the 
presence and differing material properties of the patches. Moreover, while ‘handbook’ 
values can be determined for p, D and // for a plate which is devoid of patches, those values 
usually cannot be used in the final system model with any accuracy due to nonuniformities 
in the plate and boundary conditions, variations in materials, and contributions due to 
the presence of the patches. In applications, such as that in this paper, these parameters 
are estimated using fit-to-data. techniques (see [2, 3] for the parameter estimation results 
pertaining to the plate used here). 

3. The input parameters Kf are discontinuous since they are nonzero only over the regions 
of the patches. While expressions for these constants are derived in [11], they too must 
be estimated in applications due to manufacturing variations in the patches. 

It is readily noted in the strong form of the modeling equations that the discontinuous plate 
parameters and patch input terms are differentiated, thus leading to derivatives of the Dirac 8 
‘function’. The difficulties associated with this formulation are avoided in the weak form of the 
modeling equations which is presented in the next section. 


2.2 Weak Form of Plate Model 

To provide a framework which facilitates analysis, approximation and implementation, it is 
advantageous to consider a weak form of the modeling equations. Such a formulation can be 
determined directly from Hamiltonian (energy) principles and is equivalent to that obtained l)} 7 
integration 1 >v parts after multiplication of the strong form by suitably smooth test functions. 

For the plate problem under consideration, a suitable space V of test functions is the subset 
of the Sobolev space // 2 (Fo) which satisfies the essential boundary conditions w — ~ = 0 at 
r = a. As detailed in [9, 29], a weak or variational form of the equations of motion for the plate 
is 




d 2 ij di} 

r d^od ~ m 



L > 


dll d 2 n I , 

r lh- + W‘ ,h 


fcf u i[t)\i{r , 0)X 72 i](l~f + [ gijd 
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( 2 . 2 ) 


for all i] 6 V. The overbar here denotes complex conjugation and the differential is ch/ = rdOdr. 
It is easily noted that in this form, derivatives that were originally applied to moments have been 
transferred onto test functions. This eliminates the difficulties associated with differentiating 
the piecewise constant parameters />, Z), u and Co found in the internal moments as well as the 
discontinuous input parameters Kf , • • • , Kf . 


2.3 State Approximation 

As discussed in [9, 29], an appropriate choice for the basis and Fourier-Galerkin expansion 
of the plate displacement, when considering clamped boundary conditions, is Bff{r,0) = 
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Here B™(r) is the n th modified cubic spline satisfying B™ (a) = dB %} a ^ — 0 with the condition 
= 0 being enforced when m — 0 (this latter condition guarantees differentiability at the 
origin and implies that 


N m 


N , m = 0 

N + 1 , in ^ 0 


where N denotes the number of modified cubic splines). The total number of plate basis 
functions is J\f = (2 M + 1 ){N + 1) — 1. As discussed in the [9, 29], the inclusion of the weighting 
term with 

f 0 , m = 0 

m — { 

\ 1 , rn ^ 0 

is motivated by the asymptotic behavior of the Bessel functions, which make up the analytic 
plate solution, as r — > 0. It also serves to ensure the uniqueness of the solution at the origin. 
The Fourier coefficient in the weight is truncated to control the conditioning of the mass and 
stiffness matrices (see the examples in [9]). 

To obtain a matrix system (again, see [9, 29] for a careful derivation with complete details), 
the N dimensional approximating subspace is taken to be H^ r = span{B^} and the product 
space for the usual corresponding first-order vector system is 7 ~L xr x . The restriction of the 
first-order form for the infinite-dimensional system (2.2) to the space 7 x TL^ then yields the 
matrix equation 
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where i) A ( t ) = [wf(t), wfit), • ■ • , wj(f(t)] T denotes the column J\f vector containing the approx- 
imate state coefficients (see (2.3)). The component matrices and vectors are given by 

A'n" = R’di + KD 2 + &D 3 + A'd 4 + R’d5 


— A co i + R Cd 2 + A Cd 3 + R Cd 4 + A c d 5 + f 7 Bk Bf ( h 
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The index ranges here are fc,( = 1, • • • , .V. The matrices /\ CD i — A' Cd . 5 are defined similarly with 
the inclusion of the parameter cp in the various integrals. Finally, we remind the reader that 
fK D, v and cp are piecewise constant in these definitions due to the presence of the patches. 

For application purpose, it useful to note that the matrix system for the plate can thus be 
written as the Cauchy system 

2 f(t) = A*y"(t) + B xr u(t) + </{t) 

(2 • 5, 

where y A (t) = [i) Ar (t),d Af (t)] T = [w?{t), ■ ■ - , ibf(t). ■ ■ • , wj[r{t)] T denotes the column 

2 Af vector containing the generalized Fourier coefficients for the approximate displacement and 
velocity. In this form, the control problem can be readily discussed. 


3 Continuous Time Control Problem 

In the preceding discussion leading from the infinite dimensional PDF model to the finite dimen- 
sional matrix approximation, the superscript J\f was used to denote the level of discretization; 
i.e., the number of Fourier/spline basis elements used to approximate the state. This notation 
is standard in the theory of finite element and spline approximations of infinite dimensional 
systems. In finite dimensional control theory, however, the level of discretization is typically 
fixed and these subscript JV"’s are usually omitted to simplify notation. We will do the same in 
this and subsequent sections so as to remain consistent with standard control notation. 


3.1 Initial Displacement and Velocity — LQG Control Law 

We consider first the V-dimensional systems 

y{t) = Ay(t) + Bu{t) , 2/(0) = 2/0 

2 lob(t) = Cy(t) (3.1) 

z(t) = Hy(t) + Gu(t) 



where y 0 t, denotes observations in 1R P and C is a P x Af observation matrix whose structure 
is determined by the manner and number of observations being used (the specific C matrix 
used in the plate experiments is described in Section 4.2). Moreover, ~ £ lR r denotes the 
performance output obtained under the assumption that G and H are time invariant matrices 
satisfying H T G = 0. In the event that P = M and C is an identity, the optimal control u can be 
obtained from standard linear quadratic regulator (LQR) optimal control theory. The number 
of observations P is usually limited, however, and we concentrate instead on the case P < A" 
which occurs when the full state is unavailable and must be reconstructed using a compensator 
(e.g., see [25]). 

The general control problem for this case consists of determining the voltage u which mini- 
mizes the performance index (or cost functional) 

jm = r\z(t)\\u 

JO 

= / {(Qy{i)iy{t)) + (Ru(t),u(t))} dt 

Jo 

subject to (3.1). The J\f X Al‘ matrix Q can be chosen to satisfy various design criteria including 
frequency windowing, the weighting of various state components, or minimization of certain 
energy measures. The Q matrix used here was chosen using energy considerations and con- 
struction details can be found in Section 4.2. The s x s matrix R weights the voltages to the 
various patches or patch pairs. 

Because full state information is not available in most applications, the state must be esti- 
mated or reconstructed from observations before a controlling voltage can be determined. We 
consider here a full order compensator or observer of Luenberger type [25] and refer the reader 
to [16, 22] for details on reduced order observers. 

The compensator or reconstructed state satisfies the matrix system 

i)c(t) = Ay c (t) + Bu(t) + F[y ob {t) - Cy c (t )] 

!Jc(0) = y c o 


with the optimal voltage 


u(t) = —I\y c (t) 


where F and I\ denote the compensator and feedback gains, respectively. We note that I\ 
and F are chosen so that the reconstruction error | y(t) — y c (t)\ — > 0 as t — » oo. Under usual 
observability and controllability hypotheses (see [25]), the optimal feedback and compensator 
gains are given by 

K = R~ 1 B t U 

t - i (3.2) 

F = PC t R ~ 1 ’ 


where 11 and P are unique nonnegative-definite solutions to the following feedback (regulator) 
and compensator (observer) algebraic Riccati equations 


n a + a t u - ubr~ 1 b t u + q = o 

PA t + AP- PC t R~ 1 CP + Q = 0 , 
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respectively. As was the case with the matrices Q and R, the matrices Q and R are design 
criteria for the specific control application under consideration (specific choices used in the plate 
experiments are summarized in Section 4.2). We point out that in terms of the component 
matrices and control voltage, the compensator can be expressed as 


Vc{t) = 


A - BR~ 1 B t U 


y c {t) + PC T R- x C[y(t)-y c (t )] . 


The control law just described must be implemented in real time in order to be a viable 
method for reducing vibrations in physical structures. To facilitate implementation, it is pru- 
dent to calculate offline as many components as possible and then treat those pre-calculated 
components as filters when performing online computations. The method for continuous time 
implementation is summarized and offline and online components are categorized in the follow- 
ing algorithm. 


Algorithm 3.1.1: Continuous Time Control of Initial Displacement 

and Velocity 


Offline (?) 

(*'*') 

(***) 

Construct matrices A, B, C, Q , P, Q , P 
Solve Riccati equations (3.3) for If and P 
Construct I\ = R~ 1 B T Yl and F = PC T R~ l 

Online (i) 

Collect data y 0 b(F = Cy(t) 

(«) 

Solve the ODE system 


ilc(t) = [A- BK - FC] y c (t) + Fy ob (t) 


= A c y c (t) + Fy 0 b(t) 

(Hi) 

Calculate the voltage u(t) = —I\y c (t) 


We note that the expensive (time consuming) calculation of the component matrices A, B, Q, R, 
C,Q,R and Riccati solutions II, P is performed offline with the results loaded into the control 
code as datafiles. This leaves the integration of the system y c (t) = A c y c (t) + Fy ob (t) as the 
primary computation to be performed during implementation. Issues regarding the numerical 
integration of the system as well as the effects of discrete time calculations will be discussed in 
Section 4, and a discrete version of Algorithm 3.1.1 is summarized in Algorithm 4.1.1. Finally, 
the observed system, performance index and state estimator and LQG control law for the system 
are summarized in Algorithm 3.4.1 in Section 3.4. 

While this compensator does provide the desired performance, it may lack robustness in 
some applications. In cases where added robustness with regard to certain types of system or 
observation noise and modeling errors is required, an H ™/ MinMax compensator of the type 
described in the next section can be used. 
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3.2 Initial Displacement and Velocity - H °°/ MinMax Control Law 

We consider now the design of a. dynamic compensator which is robust with respect to certain 
types of state and measurement uncertainties or disturbances (see [15] for details). To incorpo- 
rate such uncertainties, we let io(t) G 1R 9 denote input and output disturbances. The system, 
with no exogenous force, is then given by 

y(t) = Ay(t) + Bu(t) + Dw(t) 

Vob{t ) = Cy(t) + Ew{t) 

z(t) = Hy{t) + Gu(t) . 

For this discussion, we assume that the input and output disturbance are independent and 
hence DE T = 0. This is a matter of convenience and the dependent case can be handled 
similarly after slight modifications are made (see [15]). 

In this case, the min/max optimization problem leading to the controller consists of finding 
a controller u* G U = T 2 (0, oo; IFF) and disturbance w* G H = L 2 (0, oo; 1R 9 ) such that 

J* = inf sup ) = J~f(u*,w*) 

uei' 

for the disturbance augmented functional 

J 7 (u,u>) = {|o(f)| 2 -7 2 M0| 2 }^ 

= J o {(Qy(t),y(t)) + (Ru(t),u{t)) - 7 2 (w(t), w{t))}dt . 

As noted in [12, 15], the results from this optimization problem yield a bound 7 for the H^’ 
norm of the transfer function from disturbance C(w) to the performance output C{z) where C 
denotes the Laplace transform. 

Under the assumption that the pair ( A,B ) is stabilizable, (A, C) is detectable, (/l, G) is 
controllable, and (/l, H) is observable, one can prove the existence of (minimal) positive definite 
solutions II and P to the algebraic Riccati equations 

ua + a t ii - n br~ 1 b t - 7 ~ 2 g ] n + q = 0 
pa t + ap-p [c t r- 1 c - 7 ~ 2 q\ p + q = 0 

for a given attenuation 7 > 0. Moreover, if the spectral radius p of PII satisfies the condition 

p(Pn)< 7 2 , or n-7 2 P" 1 <0, 

then there exists a unique optimal controller 

«*(i) = -R~ 1 B T Uy c {t) . 

The state estimator y c {t) G 1R A( " satisfies 

y c {t) = A c y c (t) + Fy 0 b(i) 
i/c( 0 ) = y e 0 
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where 


A e = /I - BI< - FC + 7' 2 0n 


F = 


/ - 7 - 2 pn 


-1 


PC R 


Tf >- 1 


The implementation issues concerning the method are similar to those discussed in the 
algorithm for the LQG controller but determination of the Riccati solutions II, P and filter F 
are complicated by the fact that a suitable design value of 7 must be determined before matrix 
calculations can proceed. Fortunately, these calculations can be performed offline and resulting 
matrices input as data files for the online computations. Hence the actual online controller can 
run at the same rate as that obtained using the LQG methodology. 

The form of the observed system, the extended performance index and the //°°/MinMax 
control law are summarized in Algorithm 3.4.2 in Section 3.4. The various components can be 
compared with the LQG components summarized previously as well as those which arise when 
a periodic exogenous force is considered. 


3.3 Periodic Primary Excitation 

For the case in which a periodic exogenous force drives the system, knowledge of that force 
can be used extend previously discussed results to include the effects of periodicity. The Af- 
dimensional system in this case is 


y{t) = Ay(t) + Bu(t) + g(t) , 2/(0) = ij(t ) 

(3.4) 

Vob(t) = Cy(t) 

where g(t) € 1R A is periodic with period r. That periodicity is then reflected in the performance 
index 

J{u)= [ {{Qy(t),y(t)) + (Ru(t),ii(t))}dt 

JO 

which is minimized subject to (3.4). 

With A and F defined in (3.2), the reconstructed state in this case satisfies the system 
y c {t) = Ay e (t) + Bu(t) + F[y ob (t) - Cy c (t)\ + g(t) 

yc{0) = y c {T) 

with the optimal voltage given by 

u(t) = —Ky c (t) + R~ 1 B T r(t) . (3.6) 


Here r is a tracking variable defined by the system 

r(t) = -[A-BK] T r(t) + Ug(t) 
r(0) = r(r) 


where n solves the first of the algebraic Riccati equations (3.3). We point out that in this 
case, the voltage contains two contributions. The first incorporates transient information by 
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feeding back state estimates while the tracking component incorporates information regarding 
the periodic force. 

Combining (3.5) and (3.6) yields the single expression 

y c (t) = [A- BK } y c (t) + FC [y(t) - y c (t)} + BR~ l B T r(t) + g(t) 

yd 0) = y c (r ) . (3 ' 8) 

for the state estimator in terms of the tracking variable. The state estimate at time t is then 
obtained by integrating (3.7) and (3.8) after solving the necessary Riccati equations. 

As was the case when considering control of the unforced system with initial displacement 
and velocity, the computations can be categorized with respect to those which can be performed 
offline and those which must be done online. The following algorithm summarizes the continuous 
time control method for periodic excitation and categorizes the offline and online components. 


Algorithm 3.3.1: Continuous Time Control with Periodic 

Exogenous Force 


Offline (i) 

(ii) 

(Hi) 

Construct matrices A, B , C, Q , R, Q< R 
Solve Riccati equations (3.3) for If and P 
Construct K = R~ l B T n and F = PG T R- 1 

Online (i) 

Collect data y a b(t) = Cy(t) and force measurements g(t) 

(«') 

Solve the ODE systems 


Hi) = -[A-BK] T r(t) + Ug(t) 


= A tr v{t) + g(t) 


l )c(t) = (A — BK — FC]y c (t) + Fy ob (t) + BR~ l B T r(t) + g(t) 


= A c y c (t) + Fy a b(t) + A T r(t ) + g(t) 

(Hi) 

Calculate the voltage u(t) = —Ky c (t) + i? -1 B T r(t) 


In this case, the tracking equation r(t) = A tr r(t ) + g(t) must be solved before the state can 
be estimated by integrating the system y c (t) = A c y c {t ) + Fy 0 b(t) + A r i'(t) + g{t). Details 
regarding the solution of these two systems with discrete data and force measurements are 
given in Algorithm 4.1.2 in the next section. 

Finally, the governing equations and component matrices for the LQG controller are sum- 
marized in Algorithm 3.4.3 of Section 3.4. Expressions similar to those in Section 3.2 arise 
when an i/°°/MinMax compensator is considered for the problem and the resulting control law 
is summarized in Algorithm 3.4.4. We also refer the reader to [4] where details concerning the 
design of an Tf^/MinMax compensator for a structural acoustic system that is subjected to a 
periodic exogenous force is considered. 
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3.4 Summary of Systems, Performance Indices and Control Laws 

The modeling state equations, observed system, performance index and compensator design 
for the previously discussed cases are summarized in this section. Specifically, the LQG and 
'i/°°/MinMax control laws are summarized for the system with no primary force as well as 
the system subjected to a periodic exogenous force. Details regarding the construction of the 
individual components can be found in the preceding section or in the implementation section 
which follows. 


Algorithm 3.4.1: LQG Control Law for Unforced System 


(1) Observed System 

y(t) = A y{t) + Bu(t) 
y 0 b(t) = Cy(t) 
z{t) = Ily(t) + Gu(1) 

y G E Ar , u G IR S 
yob G E p 

:er, ii t g = o 

(2) Performance Index 


j(u)= rotten 

Jo 


= f {(Qy(t),y(t)) + {Ru(t),u(t))}dt 
Jo 

Q = H t H , R = G t G 

(3) State Estimator and Control Law 

y c (t) = A c y c (t) + Fy ob (t) 

A c = A- FC - BK 


k = r~ 1 b t n 
F = PC T il~ A 


JIA + A t U - JlBR- 1 B T n + Q = 0 
PA t + AP- PC T R~ l CP + 0 = 0 


Controlling Voltage 

n(t) = -Ktjc(t) 
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Algorithm 3.4.2: iJ°°/MinMax Control Law for Unforced System 


(1) 

Observed Syst 

em 








y{t) = 

= Ay{t) + Bu{t) + Dw(t) 




y 

G 1R A \ a G 1R S , 

w G IR ? 


Vob(t) 

— Cy(t) + 

Ew(t) 




J/o6 G IR ? , 

= 0 


=(t) = 

Hy(t) + Gu{t) 




- 

G lR r , II 1 ’ G = 

= 0 








3 

= H t II > 0 , 

R = G t G > 0 








3 

= dd t > 0 , 

R = EE t > 0 

(2) 

Performance Index 








J 7 (^, 

w) = r { 

N0| 2 -7> 

(*)l 2 

1 (It 







-D 

(Qy(t),y{t)) + (P 

u(t), 

«(<)) - 

7 2 <u 

>(/),«>(/))} (It 


(3) 

State Estimator and Control Law 






iic(t) 

= A e y c (t) + Fy ob {t) 








A c = 

A - FC- - 

bk + 7 - 2 (?n 








I\ = 

R~ l B T U 









F = 

I - j~ 2 PU 

_1 PC t R - 1 








UA + 

A t U-U 

BR- l B T - 7' 

’ 2 3 

n + 






PA t 

+ AP-P 

’c T ^" 1 C - 7 

" 2 3 

p + 

3 = o 





Controlling Voltage 








u{t) = 

= -Pyc(t) 








Note: y{PB) < 7 2 ( 

3i- n-o 2 P _1 

< 0 
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Algorithm 3.4.3: LQG Control Law for Periodic Exogenous Force 


(1) Observed System 

y{t) = Ay{t) + Bu(t) + g(t) 

2/(0) = 2 /(t) 

Vob(i) = Cy(t) 
c(2) = Hy(t) + Gu(t) 

(2) Performance Index 

J(u)= [ T \z(t)\ 2 dt 

Jo 

= / {{Qy(t),y(t)) + {Ru(t),u(t))}dt 

J 0 

(3) State Estimator and Control Law 

2 )c(t) = A e y e (t ) + Fy ob (t) + A r r{t) + g(t) 

A C = A- FC - BI\ 

A r = BR~ x B t 

K = R- l B T U 

f = pc T ir l 

n a + a t n - ubr- x b t u + q = o 
pa t + ap- pc t r~ 1 cp + q = o 

Tracking Equation 

r(t.) = A tr r{t) + g(t) 

Atr = - [/I - BK] T 
g{t) = Tlg(t) 


y e E Ar , u £ 1R S 


Job £ E p 

:£E r , H t G = 0 


Q = H r H , R = G t G 


Control Voltage 

u(t) = —Ky c (t) + R~ 1 B T r(t) 




Algorithm 3.4.4: H °° / MinMax Control Law for Periodic Exogenous 

Force 

(1) Observed System 

2 y(t) = A y/ ( f ) + Bu(t) + Dw(t) + y(t) y E IR-^ e 1R S , te E 1R 7 

2/(0) = ?/(r) 

Vob(t) = Cy(t) + Eiv(t) y ob G 1R P , £>£ r = 0 

z(t) = Hy(t) + Gu{t) z G \R r , II T G = 0 

Q = H t H > 0 , R = G t G > 0 
Q = DD t > 0 , R = EE t > 0 

(2) Performance Index 

= J q {m\ 2 -l 2 Ht)\ 2 }dt 

= J o {{Qy[t),y(t)) + (Ru(t),u(t)) - 7 2 {iu(t), w(t))}dt 

(3) State Estimator and Control Law 

I 'lc(t) = A c y c (t) + Fy ob (t) + A r r(t) + g{t) 

A c = A- FC - BK + ~f~ 2 Qn 
A r = BR~ x B t -~r 2 Q 

k = r~ x b t n 

F = [i - -r 2 P\[\~ X PC T R~ l 

n a + 4 T n - n br~ 1 b t - ~r 2 o\ n + q = o 
pa t + ap-p [c t r~ x c - 7 ~ 2 q\ p + q = o 

Tracking Equation 

r(t) = A tr r(t) + g(t) 

A tr = -[A-BK-r 2 Qn] T 

y(t) = Ug(t) 

Control Voltage 

u(t) = —Ky c (t) + R~ 1 B T r(t) 


Note: p{PU) < 7 2 or n - ^P~ x < 0 




4 Discrete Time Control Problem 

The control laws summarized in the last section were derived under the assumption of continuous 
time sampling of observed state and force data. Hence it was assumed that y 0 b(t) and g(t) were 
available for all i within the temporal interval of interest. Moreover, it was assumed that 
y[t),y c {t) and r(t) could be obtained through exact integration of the state, state estimator 
and tracking equations. When implementing the method, however, one has available only 
discrete data values and the differential estimator and tracking equations must be numerically 
approximated at discrete time values. The manner in which the control laws are implemented 
in discrete time and the influence of this discretization on the overall performance are discussed 
in this section. 

In the last section, it was demonstrated that in the case with no primary exogenous force, 
the state equation, state estimator and controlling voltage for the LQC1 and H™ / MinMax 
problems had the form summarized in column 1 of Table 4.1. The corresponding quantities for 
a system subjected to a periodic exogenous force are summarized in the second column. Details 
regarding the component matrices can be found in the previous section. 



No Primary Input 

Periodic Exogenous Force 

State 

Equation 

y{t) = Ay(t) + Bu{t) , 2/(0) = ijo 
Debit) = Cy(t) 

y(t) = Ay(t) + Bn(t) + y(t) , y( 0) = y(r) 
Vob(t) = Cy(t) 

State 

Estimator 

Vc(t) = A c y c (t.) + Fy ob {t ) 
Vc( 0) = y C0 

y c (t) = A c y c (t) + Fy oh (t) + A r r(t) + y(t) 
2 / c ( 0 ) = y c (r) 

Tracking 

Equation 


r(t) = A tr r(t) + y[t) 
r(0) = r(r) 

Control 

Voltage 

= -R~ 1 B T Ey c (t) 

u(t) = -R- l B T By c (t) + R~ 1 B T r(t) 


Table 4.1. State, state estimator and controlling voltage for the systems with no primary 
input and a periodic exogenous force. 


The relationship of these components to the experimental plate setup, driven by a periodic 
exogenous force, is illustrated in Figure 4.1. For that setup, the physical process consists of the 
clamped circular plate with attached piezoceramic patches. The state equations represent the 
spatial discretization of a PDE model of the process with physical parameters estimated using 
experimental data so that the model accurately captures the plate dynamics. 
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Physical Process 



Figure 4.1. Experimental plate setup, modeling equations and computations necessary for 
determination of controlling voltage. 

We point out that when modeling the process or performing simulations, the observations 
!job(t) are obtained from the approximate state through the relation y„h{t ) = Cy{t) where C 
is the matrix observation operator. When experimentally implementing the method, temporal 
data from the process is used to determine the observations y 0 b{t)- This often requires processing 
either through digital or analog filters or in the software. In the circular plate experiments, 
accelerometer data a(t) was integrated to yield velocity measurements v(t) to be used for 
the observations y 0 b(t) (see the final subsection of this section for details). Once y 0 b(t) has 
been obtained, the state is estimated by numerically integrating the estimator equation and a 
controlling voltage is calculated. This voltage is then fed back to the process. 

Because data, is processed in a digital manner, the observations y 0 b and force measurements 
y are obtained only at discrete times tj. The rate at which this data can be sampled is governed 
by the data acquisition system, the software being used and the number of calculations required 
between samples. In particular, the tracking component r(tj), state estimate y c (tj) and voltage 
u[tj) must be calculated before the arrival of data at time f J+1 . While details regarding these 
calculations are postponed until Section 4.1, it is noted here that the system sizes must be 
minimized in order to permit real-time approximation of the estimation equation. This is 
a major motivation for choosing appropriate, accurate approximation techniques to spatially 
discretize the modeling PDE’s. 

As indicated previously, the open and closed loop process dynamics can be simulated by 
using the state equations to model the plate. The state estimate and state are then approxi- 


18 




mated by simultaneously integrating the corresponding equations. The integration of the state 
equation can be performed using any sufficiently accurate ODE routine which is efficient for the 
system under consideration; for example, a variable order, variable stepsize method was used to 
solve the stiff system which arose when simulating plate dynamics. If simulations demonstrating 
the levels of control that can be obtained under “optimal” conditions are desired, the estimator 
equation can be integrated using the same high-order routine. On the other hand, simulations 
representing the attenuation levels that can be expected under “implementation” conditions can 
be obtained by incorporating values of y 0 b and g calculated at discrete times and approximating 
ijc(tj) using the techniques employed when implementing the method. Simulation results using 
both techniques can be found in [14]. 

4.1 Approximation of the Estimator and Tracking Equations 

In order to obtain tracking values and state estimates to be used when calculating controlling 
voltages, the solution to the tracking and state estimator equations must be numerically ap- 
proximated. If the goal is solely to perform simulations, this can be easily accomplished using 
the same ODE solver used to integrate the state equation (indeed, the state and estimator 
equations can be combined into a single system and integrated simultaneously). This is not 
practical when experimentally implementing the method, however, and one must typically per- 
form the work subject to the following criteria. The method must be sufficiently efficient so 
as to facilitate real-time implementation and sufficiently accurate so as to resolve system dy- 
namics. The systems are quite often stiff which implies that either ('/-stability or a-stability is 
important . Finally, the difficulties in storing past data make it prohibitive to use many popular 
multistep methods. 

For the experiments performed with the circular plate, the sample rate was sufficiently fast 
(and hence A/ was sufficiently small) that a modified backward Euler or trapezoidal method 
produced adequate results. We illustrate here such a modified backward Euler method. 

4.1.1 Initial Displacement and Velocity 

Considering first the compensator/estimator system with no primary exogenous force, we find 
that the modified Euler approximate to the solution at time tj + 1 is given by 

Vcj+i = {I — AtA c ) 1 y Cj (I — AtA c ) 1 Fy 0 b(tj) 

A c y Cj T - :yob(tj ) • 

The method is modified in the sense that current observation values y 0 b{tj) are used as input 
since futures values at tj+ 1 are unknown. The time step A t is dictated by the sample rate. We 
point out that the matrix A c = (I - A M c ) -1 and vector T c = (I - AM C ) -1 F can be computed 
offline and then loaded as datafiles for the online computations. Hence the implicit nature of 
the method, which is necessary to ensure stability, does not slow the implementation. The 
discrete time implementation of the method is summarized in Algorithm 4.1.1 which follows. 
The definitions of the component matrices can be found in Section 3, and Algorithm 4.1.1 can 
be compared with the corresponding continuous-time Algorithm 3.1.1 given in that section. 
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Algorithm 4.1.1: Discrete Time Control of Initial Displacement 

and Velocity 


Offline (?) 

{ii) 

(Hi) 

(iv) 

(v) 

Construct matrices A, B , C, Q , 7?, Q , R 
Solve Riccati equations (3.3) for II and P 
Construct K = R~ l B T Yi 
F = PC T R - 1 
A c = A - BI< - FC 

Choose appropriate A t (determined by sample rate) 

Construct A c — (I — AfA c ) -1 
T c = (/ - At A^ 1 F 

Online (?) 

Collect acceleration data a(tj) 

(«') 

Integrate to obtain y 0 b(tj) = v(tj) 

(???) 

Time step the discrete estimator system 


Vcj+ 1 = -AclJcj T P - :Vob(ij) 

M 

Calculate the voltage u(tj) = —Ky c . 


4.1.2 Periodic Exogenous Force 

The application of these control techniques to systems with both transient and steady state 
behavior involves the approximation of both the tracking and state estimator equations be- 
fore a control input is calculated. While a variety of techniques and strategies can be used 
to obtain approximate values of r(tj), which are then used when computing v/(/ ; ). these cal- 
culations must ultimately be performed in real time when implementing the method. In the 
experiments involving the circular plate, the exogenous force was measured for several periods 
and the solutions to the tracking equation were approximated and stored over a time period 
commensurate with the driving frequency. These stored tracking values were then used as a 
filter when approximating the estimated state during the remainder of the experiment. 

Illustrating with the backward Euler discretization, the approximate to the tracking solution 
was determined from the difference equation 

n + 1 = (/ - Af A tr ) _1 rj + (I - A t/Urf 1 g(tj ) 

= Atr^’j T •AtrO{tj) 

subject to the final condition ?*( T ) = 0 (when implementing the method, one can simply search 
for a “suitable” zero crossing to start the approximation). 

This approximation was continued throughout several periods of the driving force with rj 
being stored in a circular buffer. This buffer was then treated as a filter when estimating the 
state using the difference equations 

y Cj+1 = (I ~ A/A c ) _1 y C j + (I - A/T c ) -1 Fy 0 b(tj) + (7 - AfA c ) -1 A r rj + (I - AtA c )~ l g(tj) 

= Acyc } T F :yob(tj ) T A r rj + A c g{ tj ) • 
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As indicated in the following algorithm, the time intensive calculations involving matrix con- 
struction, solution of the Riccati equations and matrix inversion were performed offline prior 
to the experiments, and the matrices A tr , A c , A r and vector T c were simply loaded as datafiles. 
This, combined with the solution of the tracking filter before state estimation, yielded an al- 
gorithm which was sufficiently fast for implementation. Current efforts are aimed toward the 
simultaneous approximation of the tracking and state estimator equations during implementa- 
tion. 


Algorithm 4.1.2: Discrete Time Control of Periodic Excitation 


Offline (i) 

(*'*') 

(Hi) 


(iv) 


(v) 


Online (i) 

(«*) 

(Hi) 


(Hi) 


(iv) 


Construct matrices A, B, C, Q,R,Q,R 
Solve Riccati equations (3.3) for fl and P 
Construct K = R~ l B T n 
F = PC T R~' 

A c = A- BK - FC 
A (r = - [A - BK] T 

Choose appropriate At (determined by sample rate) 
Construct A c = (I — A tA c )~ l 

Ar = (I- A Me )" 1 A r 
Atr = (I ~ AtAt r ) 

F C = (I- AtA e )- l F 

Collect acceleration data a(tj) 

Integrate to obtain y 0 b(tj) = v(tj) 

Approximate and store tracking values 

rj+i = AtrVj + A tr g(tj) 

r(r) = 0 

Time step the discrete estimator system 

Uc H i = A c y Cj + T c y 0 b(tj) + A r rj + A c g(tj) 

y c (kr) = 0 

Calculate the voltage u(tj ) = —I\y Cj + R~ 1 B T rj 


4.1.3 Higher-Order Approximations 

In the previous discussion, a modified backward Euler method was used to discretize the state 
estimator and tracking equations. As indicated, by numerical simulations reported in [14] 
and experimental results in the next section, for small A t, this provides sufficient accuracy to 
calculate an effective feedback voltage. If more accuracy is needed, a trapezoid rule or hybrid 
method of the nature discussed on page 225 of [26] can be used. These provide increased 
accuracy without adding complexity during implementation since the components At r , A c , A r 
and T c can still be computed offline. 
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4.2 Example 1 - Matrix Construction 

The control discussion thus far has been general in the sense that it holds for general systems 
of the form 


y(t) = Ay(t) + Bu{t) y(t) = Ay(t) + Bu(i) + g{t) 

or 

Vob(t) = Cy(t) y ob (t) = Cy(t) 

as long as the pair (A,B) is stabilizable and (4, C) is detectable. Moreover, the cost functional 
matrices Q,R and observation matrices Q, R have been treated as general design criteria to be 
specified according to the application under consideration. In this example, we illustrate explic- 
itly the matrices and filters used when implementing these control techniques for a vibrating 
circular plate. 

We first note that N = 16 modified cubic splines (see (2.3)) were sufficient for resolving the 
plate dynamics in the frequency range under consideration. Due to the axisymmetric. excitation 
and response of the plate, the Fourier limit M = 0 was used in all calculations. Hence a total 
of J\f =16 basis functions were used which led to 32 coefficients in the vector y. 

The formulation and sizes of all components in the control system for the circular plate are 
summarized in Figure 4.2 and 4.3 below. 


Component 

Size 

Comments 

A = 

0 

- (a I xr ) 

i 


32 x 32 

i 

The elements comprising the 
16 X 16 matrices and 

I\f are summarized in (2.4) 

B = 

0 

(a/' v ) -1 B xr 


32 x 1 

See (2.4) for the description of 
B xr 

g(t) = 

1 1 

Jr- o 


32 x 1 

The elements of y x ( t ) are de- 
tailed in (2.4) 

C = 

1 1 

o © 

o o 

B xr (r p ,6 p ),---,B X r(r p ,e p ) _ 


p X 32 

In the experiments, accelera- 
tion data was integrated to 
obtain velocity values which 
are the second state values in 
the second-order formulation. 
Since one accelerometer was 
used, p — l. See Example 2 
of Section 4.3 for a discussion 
regarding the duality between 
control and observation. 


Figure 4.2. System and control matrices used in the circular plate experiments. 
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Component 

Size 

Comments 

Q = 

1 1 

t—* 

o ^ 

to 

O 


( 1 

> 

o Sh 

o 

1 1 


32 x 32 

A weighted mass matrix was 
used for the penalty term Q . 
As discussed in [5], this pro- 
vides a means of weighting the 
kinetic and potential energy of 
the plate. 

Example 5.1: d\ = d 2 = 1 
Example 5.2: d\ ~ d 2 = 5 

R = 

’ Ru 

Rss 


S X 8 

In the plate experiments, one 
controlling patch was used so 
s = 1. 

Example 5.1: Rn = 10” 7 
Example 5.2: Rn = 10~ 10 

Q = 

1 1 

o ^ 

rj, 

to 

^ O 

> 

! ^ 

1 1 


32 x 32 

For the plate experiments, Q 
was treated solely as a design 
parameter as compared with 
the choice in [4] where physical 
arguments were used to con- 
struct the matrix. The iden- 
tity weights were taken to be 
C!=C 2 = 1 in the experiments. 

R = 

‘ R n 

l Kpp 


p X p 

In the experiments, p= 1 since 
one accelerometer was used for 
data collection. The weight 
was taken to be Rn = 1. 


Figure 4.3. Control and observation matrices used in the circular plate experiments. The ex- 
perimental results in Example 5.1 demonstrate control of transient dynamics while the periodic 
case is illustrated in Example 5.2. 

These matrices and vectors are then employed in Algorithms 4.1.1 or 4.1.2 to create the imple- 
mentation matrices and filters which were ultimately used in the experiments. 
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4.3 Example 2 — Duality Between Control and Observation 

The duality between control and observation can be noted by considering the form of the control 
matrix B and observation matrix C. Illustrating with the case in which .s patch pairs are used 
for control, the 2Af x s control matrix has the form 


B = 



where the Af x ^ matrix B Ar has elements 


B« 


= / K>V 2 Bfx J (r,0)<l 7 
J r 0 

= / K.jV 2 Bf <h f 

J j th patch 


(here \j(r, 0) denotes the characteristic function over the j th patch). 

When data from accelerometers located at the points (ri, 0\ ),•••, {r p , 0 p ) is integrated to 
obtain velocity values, the p x 2jV observation matrix is given by 


c = [0 1 ] 


C 0 

0 c 


where 



I Bf r 6{r J ,0 j ),h 
Jro 

BfiriJj) . 


With C thus defined, it can immediately be noted that 


Vob(t) - Cy(t) 
A r 

k=i 


denotes the physical value of the velocity at the point {rj,0j) given by the state equations at 
time t. This is an approximation, to within modeling and processing error, of the actual plate 
plate velocity v(t) which is measured in experiments. 

Similarly, multiplication of the state estimator coefficients y c by C produces an estimate of 
the velocity which is then compared against the measured plate velocity when integrating the 
state estimator equation 


Vc(t) = 


A - BR- l B T n 


y c (t) + PC T R- 1 [v(t)-Cy c (i)} . 


It should be noted that as the state estimates approach the measured plate values, the estimator 
equation approaches the state equation 


y(t) = Ay(t)-BR- l B T Uy c (t) 


which is used to model the plate dynamics. 
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4.4 Integration of Experimental Data 

For the experiments involving the control of circular plate vibrations, data consisted of accel- 
eration measurements obtained from one or more accelerometers on the plate. It was then 
necessary to approximately integrate this data to obtain velocity values so as to have a state 
variable for control calculations. An issue which turns out to be crucial when approximately 
integrating experimental data concerns the robustness of the integrator with respect to inexact 
initial conditions and DC gains or biases (added constants) in the data. The inexact initial con- 
ditions can be due to unknown system contributions, static shocks during system connections, 
et cetera. While careful calibration can alleviate some of the uncertainty in initial conditions, 
it cannot fully eliminate the problem. The problem of gains or biases due to small DC voltages 
in the system can also be minimized but never fully eliminated. Hence an integrator which is 
minimally affected by uncertain initial conditions and DC gains in the data is crucial for success 
when approximately integrating data. 

Here we consider two techniques for approximately integrating acceleration data to obtain 
velocity values in accordance with the relation 


v{t) = a(t) . 

Essentially, the idea is to replace the integration by either the first-order differential equation 


v + fh> 


RC 


or the second-order equation 


v + fin + Q 2 t> = 


RC 


(4.1) 


(4.2 


(see [21]). The design parameters H and RC are frequency and time constants, respectively, 
which are chosen so that RC = 1 and u: > 6fi, where u> is the smallest observed frequency. For 
solution, (4.2) is written as the first-order system 


(4.3) 


V 


' -n 

n ' 


V 


- a(t) * 

— 


+ 

RC 

e 


—Q 

0 


e 

0 


=>i = Az + f 

where e = — f Iv. The integration of (4.3) is subject to the initial conditions 


' 1»(0) ' 


’ Vo ' 

. e(0) . 


6 o 


(4.4) 


As will be detailed in the subsequent discussion, the first-order integrator (4.1) is robust 
with respect to inexact initial conditions but propagates DC gains in the data. In the second- 
order integrator (4.2), or system (4.3), both disturbance in initial conditions and DC gains 
in acceleration data are exponentially attenuated. Moreover, the frequency response of this 
approximator is very close to that in the original signal for u> > 6fL Hence this latter method 
provides an accurate and robust means of approximately integrating experimental data. 
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4.4.1 First-Order Approximate Integrator 

Here we examine the properties of the first-order approximate integrator (4.1). We consider 
first the case in which the exact initial condition Vo is known and the acceleration a(t) is free 
from DC gains or biases (added constants). If we let C(s) = £{t>(f )} and /l(s) = £{a(f )}, then 
Laplace transformation of the system (4.1) yields 


V(s) 


1 1 
7 Tn v ° + s + n 


i 

RC 


A(s) . 


Hence the approximate integrator (4.1) is a single pole filter. Inverse transformation then yields 


«(<) = ^'-'aisUs 


(4.5) 


as the solution to (4.1). 

Similarly, if vq denotes a perturbed initial condition and a DC gain g a is present in the data, 
then the solution is given by 


v(t) = e Qt v 0 + 


RC Jo 




«(*) + 9a] ds 


It follows immediately that. 


i,{t) = v(t) + e Qt (u 0 - Vo ~ -J^j) 


+ 


9a 

rcq ‘ 


It is first noted that the perturbations in initial conditions exponentially decay with the rate 
of decay influenced by the magnitude of the parameter Q. DC gains of the order g a /(RC Q) 
remain, however, thus leading to difficulties when such biases are present in the data. Both 
properties are numerically illustrated through examples in [14]. 

The manner through which the solution (4.5) approximates the solution to the original 
relation v(t) = a(t) can be illustrated with a simple example. Consider a(t) = 1207T cos( 1207 t<). 
The solution to (4.1) for this acceleration is 


v (i) = • Q 2 ^( 190^)2 [ncos(120rt) + 1207rsin(120w<)] 


1 120ttD 

RC ' O 2 + (120tt) 2 e 


which reduces to the solution of the original relation with ft = 0 and RC = 1. For D = 167T, 
the solution v still provides an adequate approximation to the original whereas it is a very poor 
approximation with fl = 1 20 tt. This phenomenon is illustrated in [14]. 


4.4.2 Second-Order Approximate Integrator 

The second-order approximate integrator (4.2), or equivalent system (4.3), eliminates the diffi- 
culties associated with both inexact initial conditions and DC gains or biases in the data. The 
elimination of constants in the data can heuristically be attributed to the differentiation of the 
acceleration data. This can be made rigorous by analytically solving the problem. 
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We again let V(s) = £{v(f)} and A(s) = £{a(t)} and let t’o, t>i and « 0 denote initial 
conditions. Transformation of (4.2) yields 

s 2 V (s) — svq — t’i| + Cl [sV’(.s) — r>o] + H 2 V(s) = —— [s.4(s) — no] 


RC 


from which it follows that 

s + Q 


V(s) = 


.. , (vi — « 0 / RC) , 

To H — o . ^ — r"^T + 


1 


;A(s) • 


s 2 + n s + n 2 .s 2 + ns + n 2 s 2 + n s + n 2 rc 

Inverse transformation then yields the solution 
v(t) = e' m/2 cos {V3Qt/2) v 0 + ~^= sin (\/3fl</ 2 ) v 0 + sin (\Z3flt/2) ^ 

+ W? Jo e ~ Q{t ~ S)/2 cos - s )/ 2 ) ~ ^5= sin - s)/ 2) 


A(s)ds . 


As in the discussion of (4.1), we then consider the corresponding solution with perturbed 
initial conditions v 0 ,i>i and « 0 and DC gain g a . In this case, the perturbed solution v(t) is given 

by 

v(t) + e~ Qt ' 2 
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Here both the perturbations in initial conditions and the added constants in the data expo- 
nentially decay with the rate of decay dependent on the magnitude of Q. This is illustrated in 
examples given in [14]. We reiterate that while increased values of fi lead to more rapid decay of 
perturbations and biases, the solution to the differential equation less accurately approximates 
the true velocity. This, in combination with the goal of accurately preserving signal frequencies, 
leads to the condition fl < uj/6. 


4.4.3 Numerical Approximation 


Because of the potential for problems involving DC gains with the first-order formulation (4.1), 
we concentrate primarily on the second-order filter (4.2) and the corresponding system (4.3). 
We note that similar scalar techniques can be used to approximate the solution to (4.1). In 
considering numerical techniques for integrating (4.3), emphasis was placed on using a technique 
which could easily be implemented in real time. The two methods considered here are Euler’s 
method and a backward Euler’s method. The two are summarized below. 


Euler’s Method: 

z ih+i ) = [I + dtA\z(tk) + dtf(tk ) 
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Backward Euler’s Method: 


~(f*) = [/ - dtA}~' s(tk-i) + [/ - dtA]- 1 dtf(t k . _j) 


where 
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The advantage of backward Euler’s method over Euler’s method is its stability properties 
with slightly more involved matrices being the disadvantage. Numerical examples demonstrat- 
ing both methods with a variety of exogenous forces can be found in [14]. 


5 Experimental Results 

Experimental results demonstrating both the transient and steady state capabilities of the 
control methodology are presented in this section. The circular plate used in these experiments 
had a radius of 9" (.2276 m) and a. thickness of ,05 // (.00127 ???). A pair of piezoceramic patches 
having radius .75" (.01905???) and thickness .007" (.0001778???) were bonded to the center of the 
plate. In both the transient and steady state experiments, only one patch was used for control. 
In the steady state experiments, the opposite patch was used to drive the plate while it was 
allowed to remain uncharged in the transient case. The plate was mounted in a wooden frame 
by a circular aluminum collar which provided boundary conditions which were sufficiently close 
to clamped (zero displacement and slope). 

The first step in the process was the estimation of physical parameters through fit-to-data 
techniques. As detailed in [2, 3], transient plate vibrations were excited though an impact ham- 
mer strike or the input of a. voltage spike to the patches, and acceleration data was measured. 
The parameter values summarized in Table 5.1 were obtained through a least squares mini- 
mization of the difference between the model response and the measured data. These values 
were then employed when constructing the component matrices A,B and Q (see Figure 4.2) 
used during the experimental implementation of the controller. 




Physical Parameters 

p • 1 1 

Plate 

3.170 

(kg/m 2 ) 

Plate + Pzt. 

3.216 

D 

Plate 

11.151 

(N • m) 

Plate + Pzt 

11.506 

CD 

Plate 

1.443-4 

(N ■ m * sec) 

Plate + Pzt 

2.031-4 

V 

Plate 

.326 


Plate + Pzt 

.325 

7 ( sec • N/ 777 3 ) 

17.021 

IC b (N/V) - 

Controlling Patch 

.016 

IC B (N/V) 

- Driving Patch 

.017 


Table 5.1. Physical parameters used in the experiments. 


5.1 Transient Control 

To investigate the capabilities of the method for controlling transient vibrations, decaying plate 
responses generated by an impact hammer strike were considered. In each case, the strike was 
directed to the center of the plate and hence the plate response was axisymmetric. In the first 
set of experiments, data was collected from an accelerometer located at the plate center on the 
side opposite from the hammer impact; thus P r = {r\,() x ) = (0,0) in the construction of the 
observation matrix C described in Figure 4.2. The experiments were then repeated with the 
accelerometer placed at the off-center point P 2 = (ri,0i) = (2", 0) to illustrate that collocation 
between the sensor and actuator is unnecessary in this control method (see Figure 5.1 for 
sensor, actuator and impact locations). The results obtained with the off-center accelerometer 
are presented here, and the reader is referred to [14] for a discussion of the transient control 
results obtained with observations from the centered accelerometer. 



Figure 5.1. Patch, accelerometer and impact locations for the experiments involving control 
of transient vibrations. 
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Since no exogenous force was applied to the plate, the state estimator and control law 
summarized in Algorithm 4.1.1 were used to compute the controlling voltage to the patch. The 
component matrices as well as the cost functional and observation parameters used in these 
experiments are summarized in Figures 4.2 and 4.3. 

Data acquisition and processing were performed with a PC-based Texas Instruments TMS 
320-C30 digital signal processing (DSP) board. A schematic of the amplifiers, filters, DSP 
configuration and PC algorithm is given Figure 5.2. In the experiments, the accelerometer 
voltage was initially boosted by a factor of 10 and then reduced by 5 dB before reaching the 
DSP board. The controlling voltage output from the DSP was also boosted by an amplifier 
before input to the patches. This was necessary since the maximum voltage output by the DSP 
is 2.5 V whereas 60—70 V were needed at the patch. The reader is also referred to Figure 4.1 for 
an illustration of the experimental process and to [14] for details regarding the implementation 
process. We point out that the ratio |^§, illustrated in the A/D conversion, results when the 
voltage range —2.5 to 2.5 is discretized into 4096 possible digital values. The reciprocal process 
occurs when digital values are converted to analog voltages in the D/A converter. 

The control code was written in assembler in order to attain sufficiently fast sample rates 
for resolving transient frequencies excited by the hammer impact. While the code ran at rates 
greater than 7 KHz, a sample rate of 3.5 KHz was used in the experiments. This proved to be 
sufficient for resolving the three axisymmetric modes (with frequencies of 60 Hz, 227 Hz and 
512 Hz) excited in the experiments. 

Representative plots of the plate velocity (integrated from the data) at the off-center point 
P -2 in the uncontrolled and controlled cases are given in Figure 5.3 with reduction levels at times 
t = 0.5, 1.0, 1.5 sec summarized in Table 5.2. The percentage reductions were calculated by de- 
termining the ratio between the maximum values of the controlled and uncontrolled trajectories 
through one period containing the time point of interest. As illustrated by the results in the 
table and figure, the velocity level in the controlled case has been reduced by 50% before .5 sec 
and is essentially fully attenuated by 1.5 sec. We reiterate that these results were obtained with 
data obtained from an accelerometer at (2 // ,0) and a centered actuating patch thus illustrating 
that collocation is unnecessary for this control method. 

The voltage, u(tj) = — I\y Cj , was recorded in each experiment and that yielding the control 
results reported here is plotted in Figure 5.4. It is noted that the voltage has a maximum 
magnitude of 70 V. In practice, it has been observed that the patches can be used for extended 
periods at the frequencies of interest without damage or degradation of performance if the 
voltage levels are maintained below S— 10 rms V/mil [28]. Hence control voltage levels required 
for control of the transient vibrations is well within the tolerance of the 7 mil patches used in 
the experiments. 

Finally, the force delivered by the hammer impact in the uncontrolled and controlled cases 
is plotted in Figure 5.5 so as to provide a means of testing the equity of excitation levels in 
the uncontrolled and controlled experiments (the initial velocities in the two cases can also be 
compared to determine whether the same level of energy is being delivered in each experiment). 
As indicated by the results in this latter figure, the force delivered in the two cases is nearly 
identical (a slight double hit was always present when an impact hammer was used to excite 
the axisymmetric modes). 
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Time 

Acceleration 

Velocity 

.5 sec 

09. 5 

08.2 

1 sec 

88.9 

84.7 

1.5 sec 

93.8 

97.8 


Table 5.2. The percent reductions in acceleration and velocity levels at the point P 2 = (2",0) 
when feedback control is implemented. 



Plate + Patches 


Hardware Software 


Figure 5.2. Amplifiers, DSP configuration and PC Algorithm 4.1.1 for controlling a plate 
excited by an initial impact. Component matrices are defined in Figures 4.2 and 4.3. 


31 







0.15 


0.1 


0.05 


o 

Q 

to 


o 

o 

o 

> 


- 0.05 


- 0.1 


- 0 . 1 5 [- 
0 


Time (sec) 


Figure 5.3. Uncontrolled and controlled plate vibrations at (2",0) in response to an impact 
hammer hit; (Uncontrolled), — (Controlled). 


Controlling Voltage to the Patch 



Figure 5.4. The controlling voltage for acceleration data observed at (2",0). 





Hammer Forces in the Controlled and Uncontrolled Cases 



Figure 5.5. The force delivered by the impact hammer in the uncontrolled and controlled 
cases, (Controlled Case), (Uncontrolled Case). 


5.2 Control of a Periodic Exogenous Force 

A second problem under consideration concerns the control of plate vibrations when the plate 
is driven by a periodic exogenous force. To demonstrate the control capabilities in this case, a 
periodic driving voltage was supplied to one centered patch on the plate and the patch on the 
opposite side of the plate was used as the control actuator. Experimental tests indicated that 
a 350 Hz driving voltage produced a. strong plate response and all tests were conducted with 
the exogenous voltage at that frequency. 

Two sources were used to generate this exogenous signal: namely, an external oscillator and 
the PC running the control algorithm. As reported under Case 1 below, a purely steady state 
response could be considered with the oscillator-generated exogenous force since the plate was 
driven to steady state before the control program was initiated. Both a transient and stead}’ 
state response were noted in the PC-generated signal since the input of the exogenous force to 
the plate began at the same time that the control algorithm was started. This latter means of 
excitation is considered in Case 2. 

Since the system was driven by a periodic exogenous force, the discrete time Algorithm 

4.1.2 (corresponding to continuous time Algorithm 3.3.1) was used to calculate the controlling 
voltage to the actuating patch. Again, component matrices as well as control and observation 
parameters are summarized in Figures 4.2 and 4.3. 

For the ofF-center results reported here, velocity observation values were obtained by inte- 
grating data from an off-center accelerometer located at the point P 2 — (2", 0) as depicted in 
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Figure 5.1. The second input to the algorithm consisted of measurements g{tj) of the force 
driving the plate. We point out that when implementing this control method, both phase 
and magnitude information for the driving force were recjuired as compared with many other 
methods (e.g., feedforward) which require only phase information. 

In these experiments, the tracking components rj were calculated first and stored in a cir- 
cular buffer. These values were then used when calculating the state estimates ijj and voltages 
u(tj). While the implementation in this manner facilitated running the algorithm with sample 
rates on the order of 7 KHz (again, the algorithm was coded in assembler), it limits the robust- 
ness of the method with respect to changes and variations in the driving force. Current efforts 
are directed toward simultaneous solution of the tracking and estimator difference equations. 

A crucial issue when implementing the method concerns the handling of delays and phase 
shifts produced by the filters, A/D and D/A conversions, computation of the control voltages, 
et. cetera. While the amount of delay and phase shift is frequency dependent, experiments 
indicated that at 350 Hz, 30 — 40° phase shifts were introduced by the hardware. This was suf- 
ficient to destabilize the controller if left uncompensated. In the experiments, we compensated 
by first conducting an offline, numerical ‘identification’ to determine the amount of added delay 
necessary for stabilizing the controller in the presence of phase shifts of the order introduced by 
the experimental hardware. A summary of these results can be found in [14]. This is analogous 
to the online tuning or phase locking which is necessary for ensuring stability in other control 
methods. The numerical tests, summarized in [14], indicated that the introduction of a 216° 
delay in accelerometer or exogenous force data would stabilize the systems, and this was im- 
plemented in the experiments by using a DAC delay box. The compensation for phase delays 
in this manner provides merely a. first step toward optimal implementation of the method and 
one aspect of current research efforts is directed toward online compensation in the algorithm. 

Data, acquisition and processing was again performed with a PC-based Texas Instruments 
TMS 320-C30 board. A schematic of the setup is given in Figure 5.6. 
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Plate + Patches 


Hardware 


Software 


Figure 5.6. Amplifiers, DSP configuration and PC Algorithm 4.1.2 for controlling a plate 
driven by a periodic exogenous force. Component matrices are defined in Figures 4.2 and 4.3. 


Case 1: Oscillator-Generated Driving Signal 

For the results described here, the signal to the driving patch was generated by an external 
oscillator. The plate was allowed to reach steady state and then the control program was initi- 
ated. Acceleration levels measured by the accelerometer located at P 2 = (2", 0) and integrated 
velocity values for the uncontrolled and controlled cases are plotted in Figure 5.7. As noted 
from the controlled trajectories in that figure, it takes the algorithm approximately 0.0G seconds 
to calculate and store a sufficient number of tracking components rj. During that time interval, 
no voltage is fed to the actuating patch. Once the tracking calculations are completed, state 
estimation begins and the controlling voltage is computed and fed back into the system. The 
vibration levels decay for approximately 0.3 seconds and then are maintained at levels that 
are approximately 15% of those for the for the uncontrolled case for the remainder of the time 
interval. This corresponds to a 20log(rt con /rt uneon ) « —16.5 dB reduction in acceleration levels. 

While the magnitude of the controlling voltage is dependent upon the amplitude of the 
driving signal, magnitudes less than 40 V (28.3 V rms) were required to attain the levels reported 
here. At 350 Hz, this was well within the range (56 — 70 V rms in this case) that was considered 
to be safe for the patch being used (see Example 5.1 for further discussion regarding the voltage 
levels at which the patches can be driven without damage or degradation in performance). 

The results in this experiment demonstrate the effectiveness of the algorithm for controlling 
a system that has reached steady state. Hence one is not required to start the control with a 
system at rest. 
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Figure 5.7. Uncontrolled and controlled plate acceleration and velocity at (2,0) for system 
driven by periodic exogenous force; (Uncontrolled), — (Controlled). 

Case 2: PC-Generated Driving Signal 

A second mechanism for generating the driving signal is with the PC that is used to process 
data and run the control algorithm. Acceleration and velocity plots of the uncontrolled and 
controlled plate vibrations excited in this manner are given in Figure 5 . 8 . It can be seen that in 
this case, the plate starts from rest and is still being driven through a transient stage when the 
tracking calculations are completed and control begins. At that point, the controlled trajecto- 
ries are reduced to the levels noted in the purely stead}' state case whereas the uncontrolled 
trajectories are driven to steady state. Here, an 82% (15 dB) reduction in levels is noted at time 
T = 1 second. This was obtained with a controlling voltage of magnitude 1 2 V mar ( 8.5 V„ 71S ) . 
These results demonstrate the effectiveness of the control algorithm for a system undergoing 
transient oscillations before reaching steady state in response to a periodic driving force. 


Uncontrolled and Controlled Plate Acceleration 


Uncontrolled and Controlled Plate Velocity 
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Figure 5.8. Uncontrolled and controlled plate acceleration and velocity at (2,0) for system 
driven by periodic exogenous force; (Uncontrolled), — (Controlled). 
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6 Conclusions 


In this work, the experimental implementation of a PDE-based controller was considered. While 
the motivating application involves the control of vibration levels for a circular plate through 
the excitation of surface-mounted piezoceramic patches, the general techniques described here 
will extend to a variety of applications. 

For such control techniques, the first step is the derivation of a PDE model which accurately 
describes the dynamics of the system under consideration. Physical parameters in these models 
must typically be estimated through fit-to-data techniques before control applications can be 
considered. Following a brief discussion regarding the strong and weak forms of a thin plate 
model with a discontinuous control input term (due to the piecewise constant nature of the 
piezoceramic patches), continuous-time LQG and if 00 methods for systems with no exogenous 
force or a periodic exogenous force were discussed. The discrete-time approximations necessary 
for implementing the methods with digital measurements were also presented. A crucial step 
when implementing the discrete-time controllers involves the approximate integration of data 
(e.g., accelerometer data integrated to obtain velocity state values), and first and second-order 
filters for accomplishing this were discussed. Without such filters, DC biases, that are always 
present in the data, would render the integrated values useless. 

Experimental results demonstrating the control of transient and steady state vibrations were 
then presented. One advantage of the PDE-based controllers over standard frequency response 
input/output techniques is the capability for direct control of transient responses, and this was 
demonstrated in the first set of examples. A centered hammer impact was used to excite the 
plate and integrated data from an off-center accelerometer was used to reconstruct the state. 
The results demonstrate that attenuation levels on the order of 95% reduction can be attained 
by 1.5 seconds using the PDE-based controller. 

The second example demonstrates the control of transient and steady state responses when 
the plate was driven by a periodic exogenous voltage to a secondary piezoceramic patch. These 
results demonstrate that, after accounting for hardware delays, attenuation levels on the order 
of 85% were attained when control was implemented. While implementation techniques are 
still being refined, these results demonstrate the effectiveness of the PDE-based controller for 
this system and indicate the potential of these control techniques for reducing transient and 
steady state dynamics in other structural and structural acoustic systems. 
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